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METRIC RESULTS ON THE DISCREPANCY OF SEQUENCES (ana)„>i 
MODULO ONE FOR INTEGER SEQUENCES OF POLYNOMIAL 

GROWTH 

CHRISTOPH AISTLEITNER AND GERHARD LARCHER 

Abstract. An important result of H. Weyl states that for every sequence (an)„>i of dis¬ 
tinct positive integers the sequence of fractional parts of (ano)„>i is uniformly distributed 
modulo one for almost all a. However, in general it is a very hard problem to calculate 
the precise order of convergence of the discrepancy of ({anQ;})„>i for almost all a. In 
particular it is very difficult to give sharp lower bounds for the speed of convergence. Until 
now this was only carried out for lacunary sequences (an)„>i and for some special cases 
such as the Kronecker sequence ({r<t})„>i or the sequence In the present 

paper we answer the question for a large class of sequences (an)n>i including as a special 
case all polynomials an = P {n) with P G 1^[x] of degree at least 2. 


1. Introduction 

In the present paper we will study distribution properties of sequences {xn)n>i of fh® form 
Xn = {onCt} in the unit interval, where a is a given real and (an)„>i is a given sequence of 
integers, and where {■} denotes the fractional part. In particular we are interested in the 
behavior of the discrepancy of these sequences from a metric point of view. 


The discrepancy Dn of a sequence (a^n)„>i in [0,1) is given by 


Dn = sup 

0<a<fe<l 


An ([g, b)) 
N 


{b-a) 


where An {[a,b)) = #{l<n<A|x„G [a,b)}. The sequence {xn)n>i is uniformly dis¬ 
tributed in [0,1] if and only if lim Dn = 0- It was shown by H. Weyl [17] that for every 

N^oo 

sequence (an)„>i of distinct positive integers the sequence ({an«})n>i is uniformly dis¬ 
tributed for almost all a. In [1] R.C. Baker gave a corresponding discrepancy estimate: 
Let {an)n>i be a strictly increasing sequence of positive integers. Then for almost all a for 
the discrepancy Dn of ({on«})„>i we have 


( 1 ) 


NDn = O (logA)^+") 
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for all £ > 0. It is known that this result is best possible, up to logarithmic terms, as a 
generic result covering all strictly increasing integer sequences (a„)„>i (see [5]). However, 
in most cases of interesting particular sequences {an)n>i no metric lower bounds for the 
discrepancy of ({ana})„>i are known at all, and in these cases it is totally unclear how 
close the generic upper bound in ([T]) comes to the “correct” metric order of the discrepancy. 

One case where the precise metric asymptotic order of the discrepancy of ({anD})„>i is 
known is the case when (an)„>i is a lacunary sequence, i.e., if >1 + 5 for some constant 
5 > 0. For this case, where certain independence properties of the sequence ({onD})„>i 
can be used, W. Philipp [16] proved that for almost all a 

rn\ ^ ^ r NDn 

(2) —;= < hmsup , - < Cs 

4\/2 N^oo V^iVToglogiV 

holds. Even more precise results were obtained by K. Fukuyama [10], taking into account 
the number-theoretic structure of {an)n>i- If is interesting to note that the asymptotic 
result in ([2]) is in accordance with the Chung-Smirnov law of the iterated logarithm, which 
prescribes a discrepancy of order (A^loglog for the discrepancy of a “random” se¬ 

quence in the unit interval, almost surely. 

For the case = n, that is for the Kronecker sequence, it follows from results of Khint- 
chine [12] in the metric theory of continued fractions that 

(3) NDn = O (log N (log log N) 
and 

(4) NDn = (log N log log N) 

hold for every e > 0 for almost all a. Of course the same result also holds in the case when 
an is a polynomial of degree 1 in n. 

For the case an = it follows from a result of Fiedler, Jurkat and Korner [8] that 

NDn = O (iV^(logA^)3) 

holds. 

Apart from the case of lacunary {an)n>i) fhe classical case of the pure Kronecker sequence 
and from the example ({r^Q;})„>i only for a few further examples the precise 
metric order of the discrepancy of {{cinOi})n>i is known. A very interesting special example 
was given recently in [3]. Here for the hrst time an example for (an)„>i was given where 

the metric order of NDn is strictly between N^ and N^~^. Take for (an)„>i the sequence 
of integers with an even sum of digits in base 2. Then for almost all a we have 

NDn = O {N^+^) 
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and 

NDn = 

for all £ > 0, where a is a constant of the form k = 0,404.... Interestingly, the precise 
value of the constant k is still unknown; for details see [3] and [9]. 

Until now, however, to the best of our knowledge, nothing was known on the exact metric 
order of the discrepancy of the sequences for fc > 3 or related polynomial 

sequences, apart from the general upper bound in (p). In particular, it seems that no 
non-trivial lower bound whatsoever was known in this case. 

In this paper we will show that for a large class of sequences (an)„>i the discrepancy 
of {{cinC(})n>i has au asymptotic order of roughly This class in particular contains all 
sequences a„ = P{n) with P (n) G Z [a:] of degree larger or equal 2. Consequently, together 
with the results from ([3]) and (jl]), we have now a fairly complete understanding of the met¬ 
ric discrepancy behavior of sequences generated by polynomials with integer coefficients. 

Our main result is the following: 

Theorem 1 . Let P G Z [x] be a polynomial of degree d>2 and let be an arbitrary 

sequence of pairwise different integers with |m„| < n* for some t E N and all n > n{f). 
Then for the discrepancy of the sequence ({P (m„) we have for almost all a 

NDn > 

for all e > 0 and for infinitely many N. 

As a direct consequence of Theorem [T] (by choosing = n) we obtain the following 
corollary. 

Corollary 1 . Let P G Z [x] he a polynomial of degree d >2. Then for the discrepancy 
of the sequence ({-P(n)Q;})„>^ we have for almost all a 

NDn > 

for all e > 0 and for infinitely many N. 

The same estimate for example also holds if we choose = Pn, the n—th prime, or 
= [n/5] for some fixed f3 > 1. 

Theorem [T] is a consequence of the more general 

Theorem 2. Let / : N —)■ Z &e a function with \f {n)\ < N for some t G N and all 
n > n(t). Set 

^fi^) ■■= {(x^y) eN xN \ f{x) + f{y) = n} 

and 

Af{n) := {{x,y) eN xN \ f{x) - f{y) = n} 
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and assume that we have Af{n) = O (|h,|^) as \n\ oo or Af{n) = O (|n|^) as n —)■ oo for 
all e > 0. Then for almost all a for the discrepancy of the sequence ({/(^)a})„>i we have 

NDm > 

for all e > 0 and for infinitely many N. 

Theorem [2] is a consequence of the following result which implicitly already was used in [3], 
and which shows that under certain conditions on (««)„>! the metric behavior of the 
discrepancy of ({«««})„>! is determined by the Ll-norm of the exponential sums in a^a. 

Theorem 3. Let (an)„>i be a sequence of integers such that for some t E N we have 
\o-n\ < for all n large enough. Assume there exist a number r G (0,1) and a strictly 
increasing sequence of positive integers with {B')^ < < B^ for some reals 

B', B with 1 < B' < B, such that for all e > 0 and all L > L (e) we have 


( 5 ) 



da > Bf 


Then for almost all a G [0,1) for all e > 0 for the discrepancy Djq of the sequence 
({a„a})^>^ we have 

NDn > 

for infinitely many N. 


For the proof of Theorem [3] we need an auxiliary result from metric Diophantine approx¬ 
imation which is of some interest on its own. We state this result as Theorem 01 In 
the statement of the theorem and in the sequel, P denotes the one-dimensional Lebesgue 
measure. 


Theorem 4. Let {Rl)l>q be a sequence of measurable subsets of [0,1), with P {Rl) > 
for some constant B G M"*" and such that each Rl is the disjoint union of at most A^ 
intervals for some A G M"*". Then for almost all a G [0,1) for every g > 0 there are 
infinitely many integers hi with 


/in < (1 + v)^ 


1 

F(i?L) 


and 

G Rl- 


n>l 


of 


Theorem [T] seduces to the hypothesis that a polynomial growth behavior of (a„) 

degree at least 2 always implies a metric lower bound of order for the discrepancy 

of Here by polynomial growth behavior of degree d we mean not only that 

On > cn'^ for some c > 0 and all n large enough, but that the stronger local condition 

^n+1 C 


> 1 + 
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holds for some c > 0 and all n large enoughQ However, this hypothesis does not hold as is 
shown in the following Theorem O 

Theorem 5. For every integer d > 1 there is a strictly increasing sequence (an)„>i of 
integers with 

On+l ^ 1 I ^ 

On d 

for some c > 0 and alln eN such that for almost all a the discrepancy Djsi of the sequence 
satisfies 

ND^ = 0 {{log 

for all £ > 0. 

We conclude this section with some remarks on our theorems and some open problems. As 
noted before, by the Chung-Smirnov law of the iterated logarithm for a “random” sequence 
in the unit interval the quantity ND^ typically is of order (A^ log log thus, by ([I]) 

and by Theorems d] and [H roughly speaking the sequences investigated in the present paper 
exhibit nearly random behavior. Our method only allows us to obtain discrepancy results 
with an error of order A^^; it would be very interesting to improve these error estimates to 
logarithmic terms, such as those in ([1]). 


The fact that until now lower bounds have been almost non-existent in metric discrep¬ 
ancy theory comes from the fact that such lower bounds cannot be directly deduced from 
the second Borel-Cantelli lemma, which requires independence. In contrast in the hrst 
Borel-Cantelli lemma, which is used to prove upper bounds, only estimates for the size of 
exceptional sets are required, which can be deduced from simple moment estimates and 
Markov’s inequality. The lower bounds for the discrepancy of lacunary sequences come 
from an approximation of {{anOi})n>i by an independent random system, which is possible 
because of the fast growth of (a„)n>i. This is not possible for slowly growing (a„)„>i as in 
the case of the present paper. Instead we have developed a method which uses estimates 
for the L 1-norm of exponential sum plus an quasi-independence property of the dilated 
functions ({ha„a}) for h = 1,2,.... In the present paper the required lower bounds for 
Ll-norms are deduced from upper bounds for L4-norms, which can be obtained by simply 
counting the number of solutions of certain linear equations (see Section [3]). The quasi¬ 
independence property is established using a transition from norm-estimates for sums of 
dilated functions to certain sums involving greatest common divisors (GCD sums), to¬ 
gether with recent strong bounds for such GCD sums (see Section [2]). We believe that 
the relevance of this method goes far beyond the applications in the present paper, as the 
hrst general method for proving lower bounds in metric discrepancy theory beyond the 
well-known and totally different methods for lacunary sequences. We also want to point 


^Note that for a„ = n‘^ we have 


1 + 




for all n large enough. 
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out that these topics are related to a problem in the metric theory of Diophantine approx¬ 
imation posed by LeVeque [12] which is still unsolved; see the last section of |3] for details. 


The remaining part of this paper is organized as follows. In Section [2] we hrst prove 
Theorem H] and then Theorem [3l In Section [3] we deduce Theorem [2] and Theorem [1] from 
Theorem |3l Finally in Section H] we prove Theorem |5l 

2. Proofs of Theorem 0] and Theorem [3] 

Proof of Theorem Let (a) denote the indicator function of the set Rl, extended with 
period 1, and ]Il(d) := 1l(«) — Jq 


We have 


and 


1l (a) da = IP {Rl) 


1l (a) da = 0, 


and, by assumption, for the total variation Var of II on [0,1) we have 
(6) Var < 2A^. 

Furthermore we have 


(7) 

Let Hi : = 

( 8 ) 

and 

(9) 

for L large enough. 


i||^= / {Ii{a)fda = F{RL){l-F{Ri))<F{Ri). 


12 

\L 1 


(1 + ? 7 ) • Then by dehnition 


lL{ha)= [ lL{u)du+ IlL(ha) 

h<HL h<HL h<HL 

rl 1 

V / lLiu)du = HLF{Ri)>-il + T])^ 

h<Hr, do ^ 


We write 

OO 

1l (a) ~ y^ {uj cos 2Tija + Vj sin 27ija) 
i=i 

for the Fourier series of Ii. From ([3]) and a classical inequality for the size of the Fourier 
coefficients of functions of bounded variation (see for example [IB], p. 48) we have 


\Ui 


Var li 

< - - < — 

2i j 


and 


( 10 ) 
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We split the function into an even and an odd part (that is, into a cosine- and a sine- 
series). In the sequel, we consider only the even part; the odd part can be treated in 
exactly the same way. Set Gl '■= {AB (1 -|- 17 )) , where A and B are the constants from 
the statement of the theorem. Let Pl^ol) denote the G^-th partial sum of the Fourier series 
of the even part of II, and let tl^q) denote the remainder term (for simplicity of writing 
we assume that Gl is an integer). Then by Minkowski’s inequality we have 


( 11 ) 


Hl 

h=l 


(even) 


{h- 


< 


Hl 


^PL{h-) 


h=l 


+ 


Hl 


J2^L{h- 


h=l 


Furthermore, ca, Minkowski’s inequality, and Parseval’s identity imply that 


( 12 ) 


Hl 

^rL(/i-) 

h=l 2 


< HiWrLh 


< 


< 

< 

< 


Hr 


\ 


E 

i—C l+1 


^2L 


12 


HlA^ 

gT 

(1 -I- ri)^A^ 

1 . 


To estimate the hrst term on the right-hand side of fllll) , we expand pi into a Fourier series 
and use the orthogonality of the trigonometric system. Then we obtain 


'B.pdh-) 


(13) 


Hl Gl 

E E 

fcl,fc2 = lil J2 = l 


Uj^Uj^ 


jlhi=j2h2 

Gl 

Uj^Uj^ 


E 

jl,j2 = l 


#{('11, 


^ 2 ) : 1 < hi,h2 < Hl, jih 



To estimate the size of the sum on the right-hand size of flT^ . we assume that ji and j 2 
are fixed. It turns out that we have jihi = ^ 2^2 whenever 


hi = w 


32 

gcd(ii, j2)’ 


h 2 = w 


3i 

gcd(ji,j2) 


for some positive integer w 
































CHRISTOPH AISTLEITNER AND GERHARD LARCHER 


(see Section 5 of [3] for a more detailed deduction). As a consequence we have 

l<hi,h2< Hl, jihi= j2 
= > 1 : w < min 

Hl gcd{ji,j2) 


Hl gcd(, j 2 ) Hl gcd(ji, js) 


32 


3i 


< 


max(ji,j 2 ) 

Hl gcd{ji,j2) 

VJiH 


Combining this estimate with flT^ we obtain 


(14) 


Hl 

h=l 


^ Gl 

2 ii>i 2 =i 


\uj^Uj^\ gcd(ji,^ 2 ) 
2 \/7c2 


The sum on the right-hand side of the last equation is called a GCD sum. It is well-known 
that such sums play an important role in the metric theory of Diophantine approxima¬ 
tion; the particular sum in flTT)) probably appeared for the hrst time in LeVeque’s pa¬ 
per [15] (see also J and 0). A precise upper bound for these sums has been obtained 
by Hilberdink Hilberdink’s result implies that there exists an absolute constant Cabs 

such that 

gcd(il,i2) __ ( CabsVlog(G'T) 1 V^.,,2 

2 vTT A 

Combining this estimate with ([7|) and flTT|) (and using Parseval’s identity) we have 




Hl exp 


Cabs A/log (Gl) 


< (1 + 77 )^ exp 


VloglogGi 

CabsVlog(<^i 


P(^Z 


Vlog log Gl 


and, together with (ITTll and (IT^ . and with a similar argument for the odd part of II, we 
obtain 


(15) 


Hl 


h=l 


< (1 + 77 )^ exp 


Cabs\/log(GL) 
Vlog log Gl 


^The upper bounds for the GCD sums in m are formulated in terms of the largest eigenvalues of 
certain GCD matrices; since these matrices are symmetric and positive definite, the largest eigenvalue also 
gives an upper bound for the GCD sum. This relation is explained in detail in [5]. 
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By Chebyshev’s inequality we have 

Hl 


P a e [0,1) : 




h=l 


> (log//L) 


Hl 

h=l 


< 


(logh^n) 


and since grows exponentially in L these probabilities give a convergent series 

when summing over L. Thus by the Borel-Cantelli lemma for almost all a only hnitely 
many events 


Hl 

h=l 

happen, which by ffTSjl implies that 

Hl 

h=l 


a] 


> (loghfn) 


M 




h=l 


a] 


< (1 + 77)^/^exp 


Cabs^/L 

\/logL 


for some absolute constant Cabs- Comparing this upper bound with 
conclude that 

Hl 

liiha) S> (1 + r])^ as L —)■ oo 

h=l 

for almost all a. In particular we have 


and using ([9]) we 


oo Hl 

L=1 h=l 


a) = OO 


for almost all a, and the result follows. 


□ 


Proof of Theorem 0 By the Koksma-Hlawka inequality (see for example [6] or m) for 
every positive integer H we have 


(16) 

Let 


N Dn P- 


4H 


N 

E' 

n=l 


,2'KiHana 


h (a):= 


Bl 

E- 

n=l 


, 2'7T'tG/q^ 


We will show that for any given e > 0 for almost all a there are inhnitely many L such 
that there exists a positive integer with 


(17) 


T—/l {hia) > B 
hr 


T—e 
L 5 


which together with (1T6|) implies the conclusion of the theorem. 
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Let £ > 0 with e < | be given. From the dehnition of (a) and the fact that \an\ < n* 
for large n it is easily seen that for sufficiently large L we have 

IfiM-fiMl < 2715^5^ |ai - asl 

(18) < 27r |ai — 02 ! . 

Now let gi (a) be the function dehned by 


9l (a) := /l (j 


, -L 


for a G 


i (B'+Tb(j + i) (B'+T'") 


and for j = 0,1,..., — 1 

(for simplicity of writing we assume that is an integer). Then by fflSj) we have 

\gL - hi < 27r, 

which means that it is sufficient to prove flTTD with fi replaced hj ql. By construction 
the function can be written as a sum of at most (1?^+*)^ indicator functions of disjoint 
intervals. 


Let 


and 


Q: = 


l-T + 2e 

3i 


+ 1 


Ar.= Bl-^^Bf^ for z = 0,l,...,g. 

1 < Ao < flW 


Note that Aq = B'[ and Bj, < Ao < Dehne 


Ml := jo; e [0,1) with A* < \gL (a)| < A^+i 

for z = 0,..., g — 1. Then by (jSD we have 
Q-i , , / Q-l 


Y, P (m®) A,+i + ( 1 - ^ P (Mf) 


i=0 


Ao > 


i=0 


> 


\gL {a)\da 


Br 


for sufficiently large L. 
By Ao = B 


I we have 


^p(Afl->)A.+i> 

i=0 


Br 


for L large enough. Consequently for every L large enough there is an G {0,..., g — 1} 
with 
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which implies that 




> 


*L+1 


4gs 


l—r+Ae — 


Ql—T+be 


for L large enough. Note that, as a consequence of the construction of gi, the set is 


always a union of at most (5^+*)^ intervals. By Theorem H] we conclude that for almost 
all a for all r; > 0 there are inhnitely many integers hi with < (1 + r;)'* 


1 


{hio} G M^l^\ Consequently for almost all a we have 


P M 


(T) 


and 


hj 


\gL{hLOi)\ > 


> 


> 


(Mp>) 

1 1 Br^ 


{l + gfBl^ 4Q 


for g small enough and for infinitely many L, since by assumption B^ > {B')^ for some 
constant B' > 1. This proves the theorem. □ 

3. Proofs of Theorem [2] and Theorem [T] 

Proof of Theorem [M By Theorem [3] it suffices to show that 

rl 


N 

n=l 


^2nif{n)a 


da > N 2 ' 


for all e > 0 and N large enough. 


Using a standard trick, by Holder’s inequality with 

2 

N 


F{a) := 


E' 

71=1 


^27rif{n)a 


and G (a) := 


N 

E- 

n=l 


^27rif{n)a 


and with p = | and q = 3 we have 


N 

E' 

n=l 


,27ri/(n)a 


da = F (a) G (a) da 

Jo 


< 


{F{a))^da]" {G{a)fda^ ' 
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hence 

(19) 


Now 


and 


N 

g27ri/(n)a 

n=l 


3 / pi 


da 


N 

n=l 


^2nif{n)a 


da 


N 

n=l 


,27rz/(n)Q 


da > 


N 2-Kif{n)a 


E l\ 

n=l^ 


2 \ 2 


do; 


27zif{n)c 

Z-^n=l 


1 • 

4 \ 2 


da 


N 

n=l 


,27r2/(n)Q; 


da = 1 > 


l<m,n<N, 

f{m)=f{n) 


N 

E' 

/l=l 


j27T2/(n)Q; 


da = 1 

l</c,/,m,n<7V, 

fik)+f{l)=f{m)+f(n) 


where snmmation in the pennltimate snm is over all integers a G [—2A^*,2A^*] such that 
there exist k, I with 1 < k,l < N and f{k) + /(/) = a, and summation in the last sum is 
over all a G [—2A^*, 2A^*] such that there exist k, I with 1 < A;, / < and f{k) — f{m) = a. 
By assumption we either have 

< ( max ^ { 74 ^( 0 )} j {a : 3 1 < A:, / < A^ with/(A:) +/(/) = a} 

< N'^ 

< c{e)N^^^ 

for all e > 0 and all N > N{e), or we have the corresponding estimate for 


Hence by flT^ for all e > 0 and for A^ > A^ (e:) we have 


N 


^27rif{n)a 


da > 


nI 


1 Ar--- 

--A^2 2 


n=l (c(e)A^ 2+£)2 c(g )2 

and the result follows. 

Proof of TheoremUl By Theorem [2] it suffices to show that /(n) := P {mn) satishes 


□ 
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for all e > 0. 


Let first rUn = n. By assumption / is of degree d > 2, hence there exists a non-constant 
polynomial q E Z [x,y] such that /(x) — f{y) = (x — y)q{x,y). So, if /(x) — f{y) = n for 
some non-zero integer n, it follows that x — y is a divisor t of n, hence q{x,x — t) = j. 
This last equation has at most d — 1 solutions x. Since n has O (|n|*^) divisors t for all 
£ > 0, the assertion Af{n) = O {\n\^) follows. 

For arbitrary (m„) the result follows trivially from the special result for = n. □ 


4. Proof of Theorem [5] 

For the proof of Theorem[5]we will make use of the following Lemma on continued fractions. 

Lemma 1. For x G M let bm{x) denote the m-th eontinued fraction coefficient of x. For 
any integer ft >2 and L eN let 

L L 

SL{a) 

k=l m=l 

Then for almost all a we have 

SL{a) = O {L^+^) 

for every e > 0. 

Proof. This is a special case of Lemma 2 in [T3] . □ 

Proof of Theorem 0 For given integer d > 1 we now construct a sequence (hn)„>i which 
satishes the properties stated in Theorem [5l 

Let {bn)n>i be the strictly increasing sequence of integers running through the integers 

2^dk _|_ j<2dk-\-d—k 

for d = 0,..., 2^ — 1 and h = 1, 2, 3,.... 


We have 


Jn+l 


<2dk _j_ ji^dkA-d—k _j_ i^dkA-d—k 


= 1 + 


‘2dk _|_ j2^dk-\-d—k 

1 


2k-d j 


> 1 + 


-|- j^dk+d—k'^ d 


1 

h d 

L'r). 
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since (as is easily checked) 

2’^-^ + j < 

for j = 0 and j = 2^ and hence for all j = 0,..., 2^. So (&n)„>i has polynomial growth 
behavior of degree at least d. 

Let N be given and /cq and jo snch that for some jo E {O,..., 2*^° — l} . 

By standard techniqnes from the theory of nniform distribntion theory for the discrepancy 
Dn of ({6na})„>i we have 

ko 

NDn < 

k=l 

where by we denote the discrepancy of the snbseqnence 

({2'^^a + j2^^-^^^2^a}) ; j = 0,1 ,..., 2^ - 1. 

For this seqnence it is also well known (see for example [13]) that 

2k 

< Cabs bm ((2(''-i))"2''a) , 

m=l 

with an absolute constant Cabs, hence 

2ko 2ko 

N D]\j ^ Cabs E E '>■» ('3'‘2‘'“). 

k=l m=l 

where (3 = 2^'^“^^ and b^ix) again denotes the m-th continued fraction coefficient of x. By 
Lemma 1 therefore for almost all 2'^a (and consequently also for almost all a) we have 

NDn = O {kl+^) = O ((log as N ^ oo 

for all e > 0 (note that N > 2^°). This proves the theorem. □ 
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